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ABSTRACT 


Thr««  dlm«nsional  Dlrlohlat  problemi  for  Au«F(u)  , 

>  0,  ar«  treatad  numarloally  by  an  axoaptlonally  fast, 
axoaptlonally  aoourate  numerical  method*  Programming  d»> 
tails,  numerous  examples,  and  mathematioal  theory  are  supplied. 
Extension  of  the  method  in  a  natural  way  to  n-dlmensional 
problems  is  indicated  by  means  of  a  4-dimensional  example. 


ON  THE  APPROXIMATE  SOLUTION  OF  Au  »  F(u) 


D.  Greenspan  and  M.  Yohe 


1.  Introduction.  Because  of  its  importance  in  such  fields  as  potential  theory, 
automorphic  functions,  and  electron  radiation  [Ij-Es],  this  paper  will  be  concerned 
primarily  with  Dirichlet  problems  in  three  dimensions.  It  will  be  assumed  that 
the  problems  need  not  be  reducible  to  ones  in  two  dimensions  so  that  special 
assumptions,  like  axial  symmetry,  will  be  precluded. 

Precisely,  we  will  consider: 

Problem  D.  Let  G  be  a  three  dimensional  closed,  bounded,  simply  connected 

3 

(contractible)  domain  in  E  whose  interior  is  R  and  whose  boundary  is  S  . 

For  all  real  u,  let  F(u)  be  defined  and  differentiable  and  let 

(1.1)  F  >0  . 

u  — 


Then  if  ^x,  y,z)  c  C(S),  find  a  function  u(x,  y,  z)  which  is  a  solution  on  R 
of 


(1.2) 


9^u  8^u  8^u  . 

—  +  -Y  +  — =F(u) 
8x  8y  8z 


and  which  satisfies  both 


(1-3)  u  *  (x,  y,  z)  €  S 

(1.4)  U€C(G)  . 


Sponsored  by  the  Mathematics  Research  Center,  United  States  Army,  Madison, 
Wisconsin  under  Contract  No.:  DA-11-022-ORD-2059. 
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Under  quite  general  assumptions  on  S,  such  as  possession  of  the  cone 
property  [1,  p.  233],  It  Is  known  that  Problem  D  has  a  unique  solution  [1,  p.  372], 
and  It  Is  only  with  such  problems  that  we  shall  be  concerned.  But  since  It  Is 
not  known.  In  general,  how  to  give  this  solution  in  closed  form,  attention  will 
be  directed  toward  approximating  it.  A  finite  difference,  digital  computer 
technique  which  has  proved  exceptionally  fast  and  exceptionally  accurate  will 
be  described  and  both  mathematical  and  experimental  support  for  the  method 
will  be  pro’,  xuv  . 

Extension  of  the  method  in  a  natural  way  to  any  number  of  dimensions  will 
be  indicated  by  means  of  a  four-dimensional  example  (a  two  dimensional  example 
already  having  been  given  elsewhere  [6]). 
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2.  The  Numerical  Method.  Let  (x,  y,  z)  be  an  arbitrary  point  of  G  and  let 
h  be  a  positive  constant.  The  set  of  points  (x  ^  ph,  y  4  qh,  z  4-  rh)  ; 

p  =  0,  41,  42,  . . . ;  q  s  0,  41,  42,  . . . ;  r  s  0,  41,  42,  . . . ,  is  called  a  set 

of  grid  points.  Two  grid  points  are  said  to  be  adjacent  if  their  distance  apart 
is  h  .  The  set  of  all  lines,  each  one  of  which  contains  at  least  one  pair  of 
adjacent  grid  points,  Is  called  a  lattice.  Denote  by  those  points  which 
are  either  grid  points  in  G  ^  are  points  of  intersection  of  S  and  the  lattice. 

If  (x,  y,  z)  c  [S  then  (x,  y,  z)  is  called  a  boundary  lattice  point  and 

the  set  of  all  boundary  lattice  points  is  denoted  by  .  The  set  of  all  points 
of  Gj^  which  are  not  elements  of  Sj^  is  called  the  set  of  Interior  lattice  points 
and  is  denoted  by  Rj^  . 

If  Gj^  consists  of  n  points,  we  shall  number  these  in  a  one-to-one 
fashion  with  the  positive  integers  1,  2,  . . . ,  n  and  if  (x,  y,  z)  c  G^  and  has 
been  numbered  t  ,  then  u(x,  y,  z)  will  be  denoted  by  u^  . 

Method  D.  At  each  point  (x,  y,  z)  c  Sj^,  u(x,  y,  z)  =«^(x,  y,  z)  so  that  the 

exact  solution  is  known.  At  each  point  (x,  y,  z)  c  Rj^,  let  (x+h^,  y,  z)  , 
(x-h2,y,  z),  (X,  y  +  h^,  z),  (x,  y-h^,z),  (x,  y,  z  +  h^),  (x,  y,  z-h^)  be  those 
points  of  Gj^  which  are  nearest  to  (x,  y,  z)  in  the  positive  and  negative  x,  y 
and  z  directions,  respectively  (consult  Diagram  2.1).  Of  course  0<h^^h  , 

1  =  1,  2,...,  6  .  If  then  (x,  y,  z),  (x  +  h^,  y,  z),  (x  -  h^,  y,  z),  (x,  y  +  h^,  z) 
(x,  y-h^,  z),  (x,  y,  z  +  h^,  (x,  y,  z-h^)  are  numbered  0-6,  respectively, 
then  write  down  the  difference  equation 


I£  consists  of  m  points,  appllcaUon  of  (2.1)  exactly  once  at  each  point  of 
yields  a  system  of  m  alg^alc  equaUons  In,  say,  u^,  . . . ,  ,  the 

solution  of  which,  say,  U^,  . constitutes  the  numerical  solution. 

The  final  step  then  Is  to  solve  this  algebraic  system. 


Diagram  2. 1 
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3.  Pt 

The  derivation  of  (2.1)  is  a  straightforward  generalization  of  that  given  for 
the  two  dimensional  analogue  [6]. 

In  the  important  case  when  PeO,  the  algebraic  system  generated  in  Method  D 

is  linear  and,  by  means  of  the  elementary  techniques  of  [6] -[9],  it  is  easily  shown 

that  the  algebraic  solution  vector  exists  and  is  unique,  that  linear  over-relaxation 

converges,  and  that  the  numerical  solution  converges  to  the  analytical  solution 

in  a  suitable  class  of  functions.  The  general  case,  when  F  need  not  be  identically 

zero  and  the  algebraic  system  need  not  be  linear,  can  be  supported  in  a  fashion 

completely  analogous  to  that  of  Bers  [10]  ,  But  though  existence  and  uniqueness 

of  the  solution  of  the  algebraic  system  do  follow  and  convergence  of  the  numerical 

2 

to  the  analytical  solution  for  solutions  of  class  C  (G)  can  be  established,  one 
can  rely  at  all  times  only  on  the  extended  Llebmann  method  for  solving  the 
nonlinear  algebraic  system.  As  yet,  no  general  theoretical  basis  has  been 
developed  to  support  nonlinear  overrelaxation. 
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4.  P»  In  this  section  we  give  programming  details  for 

nmnlng  Method  D  with  FaO  .  Only  minor,  natural  modifications  are  necessary 
to  extend  the  program  to  the  case  F4  0  and  to  n-dimensional  problems. 

The  program  will  require  the  following  information: 

(a)  the  mesh  size  h  .  (h  may  be  different  for  the  x,  y,  and  z 
directions  if  desired;  we  will  assume,  however,  that  the  same  h 
is  to  be  used  for  all  directions). 

(b)  the  base  point,  (x,  y ,  z ),  chosen  so  that,  for  all  points 
(X,  y,  z)€G,  x<x,  y^y,  and  z  ^z  . 

(o)  integers  Lx,  Ly,  and  Lz  such  that  the  point 

(x»  y>  2)  =  (x  +  (Lx  - 1)  h,  y  +  (Ly  - 1)  h,  z  +  (Lz  -  IJh}  saUsfles  the 
condiUon  x  ^x,  y  <,y,  z  <  z  for  all  points  (x,y,z)  «  G  . 

(d)  the  6  functions  B^,  defined  for  all  points  of  the  form 

(x  +  n  h,  y  +  n  h,  z  +  n  h),  0<n^<L^,  0<n  <L,  0<n  <L 
^  y  z  3c  X  *  y  y'  "“2  z 

If  a  point  (x,y,  z)  is  in  G,  then  Bj^(x,y,  z)  gives  the  distance 

from  (x,  y,  z)  to  S  along  the  ray  passing  through  point  i  (see 

Diagram  2.1).  If  (x,  y,  z)  4  G,  Bj(:^  y>  2)  must  be  negative  for  at 

least  one  i  . 

(e)  the  boundary  value  function  ^(x,  y,  z) 

(f)  the  7  formulas  for  the  coefficients  of  u.,  u,,  u,,  u  ,  u  .,u_,u,  in 
the  difference  equation  (2.1). 

(g)  the  over-relaxation  factor,  u  . 

(h)  the  convergence  criterion,  c  . 
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The  program  will  analyze  this  information,  calculate  the  coefficients  of  the 
difference  equations  and  set  up  an  efficient  procedure  for  the  iterative  solution 
of  the  system  of  difference  equations.  The  program  will  then  perform  the  iteration, 
and,  finally,  print  the  results. 

Before  discussing  the  particulars  of  the  program,  an  example  is  in  order. 

For  the  solution  of  Example  1,  section  5,  the  following  infcxmation  was  supplied: 

(a)  h  =  0.1 

(b)  (X,  y,  z)  =(0,0,0) 

(c)  L  =  L  =  L  =11 

'  ^  X  y  z  / 

(d)  Bj  =  sgn(0^  -  X,  where  0^  =  l-(y^  +  z^);  sgn(0^  J  +1  ° 

I  -1  if  0  <0 

B^  =  X  ^ 

B^  =  sgn  (0y)  V|0^[  -y,  where  0^  =  l-(x^  +  z^) 

j 

B4  =  y 

Bg  =  sgn  (0^)  'V|0^|  -z,  where  0^  =  l-(x^  +  y^) 

B^  =  z 

2  2 

(e)  ^pCx,  y,  z)  =  x  +  2y  -  z 

(f)  the  7  formulas  as  given  in  equation  (2. 1) 

(g)  u  =  1.  8 

(h)  €  =  10“^ 
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Although  the  method  in  no  way  depends  on  the  computer  being  used  to  solve 
the  problem  (as  long  as  enough  high  speed  storage  is  available}  or  on  the 
programming  language,  we  will  assume  for  the  sake  of  clarity  that  the  programming 
will  be  done  in  FORTRAN  II  for  the  Control  Data  Corporation  model  1604  computer. 

The  item  of  major  concern  to  us  is  the  manner  in  which  arrays  are  stored  in 
the  computer's  memory.  If  we  are  given  a  3-dimenslonal  array  A  of  dimension 
m  X  n  X  p ,  we  will  give  the  FORTRAN  program  a  DIMENSION  statement  as 
follows: 

DIMENSION  A{M,  N,  P) 


where,  of  course,  the  actual  integers  involved  must  be  specified.  If  we  want 
to  reference  a  specific  element  of  A,  say  we  use  the  following  formula: 

q  =  i  +  m  •  (j -1+ n»  (k-l))  ; 


we  reference  a^^^^  by  referring  to  the  storage  location  A  +  q  .  The  number  q  , 
calculated  as  above,  will  be  called  the  INDEX  of  the  point  a  in  the  array  A  . 

IJ  Jv 

(we  note  that  q  is  uniquely  determined  by  the  subscripts  i,  j ,  and  k)  . 

Memory  must  be  allocated  for  storage  of  the  following  data: 


(a)  two  arrays  of  dimension  L  X  L  X  L  ,  one  of  which  will  be  used  to  store 

X  y  2 

the  solution  vector,  and  the  other  of  which  will  be  used  to  control  the 


iteration  procedure. 


(b)  the  coordinate  vectors  V  ,  V  ,  and  V 

X’  y’  z 

(c)  seven  arrays,  one  for  each  coefficient  of  the  difference  equation  (2,1)  , 
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Coefflolents  will  be  stored  in  sets;  that  is,  the  set  consisting  of  the 
i^  element  from  each  of  the  seven  coefficient  vectors  will  comprise 
the  coefficients  of  one  difference  equation.  Duplicate  sets  will  be 
stored  only  once,  and  several  distinct  equations  may  well  use  the  same 
coefficient  set.  The  number  of  sets  depends  on  the  problem;  it  is 
suggested  that  the  pogram  be  given  an  upper  limit  and  that  a  check  be 
made  to  see  that  this  limit  is  not  exceeded.  It  is  advisable  to  point 
out  the  number  of  sets  actually  computed  by  the  program;  this  gives  a 
basis  for  estimation  for  future  problems  as  well  as  providing  a  check  on 
the  poblem  being  solved. 

The  pogram  starts  with  the  values  x,  y,  and  z ,  and  generates  the 

coordinate  vectors.  It  then  examines  each  point  of  the  form  (x+n  h,y-hn  h,  z-i-n  h) 

X  y  z 

to  determine  whether  it  is  in  the  region  R  .  (The  functions  are  used  to  make 
this  determination).  If  the  point  is  in  R,  the  coefficients  of  the  corresponding 
difference  equation  are  computed,  the  solution  vector  entry  is  initialized,  and 
the  appopiate  control  information  is  computed  and  stored.  If  the  point  is  not 
in  the  region  R,  the  solution  vector  entry  is  set  to  1.0,  and  the  control  entry 
is  made  negative. 

After  this  analysis  is  complete,  we  have  in  storage  the  following  Information: 

(a)  The  "solution  vector"  array  U  .  The  entries  u^^j^  are  initialized  to 
zero  if  a  difference  equation  is  to  be  solved  for  the  corresponding  point; 
othpwlse  u^jj^  mjist  contain  1.0  . 
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(b)  The  control  array  KON  .  This  array  is  the  heart  of  the  computational 
procedure.  KON  is  negative  if  no  difference  equation  is  to  be 
solved  for  the  corresponding  point.  If  a  difference  equation  is  to  be 
solved  for  the  corresponding  point,  then  KON  contains  two  items 
of  information:  the  lower  half  of  the  memory  cell  contains  an  integer 
which  specifies  the  index  i  of  the  NEXT  point  at  which  a  difference 
equation  is  to  be  solved;  the  upper  half  of  the  memory  cell  contains 
the  Integer  J  which  specifies  the  coefficient  set  to  be  used  in  solving 
the  NEXT  difference  equation.  If  the  point  in  question  is  the  last  point, 
the  entry  is  zero.  The  information  for  the  first  point  is  stored  in  KON(0} 

(c)  The  coefficient  vectors;  these  vectors  contain  the  coefficient  sets  to  be 
used  in  computing  the  solution  at  the  various  points.  This  technique 
of  keeping  only  the  unique  coefficient  sets,  suggested  by  D.  Van  Egeren, 
effects  a  considerable  saving  in  high  speed  storage  requirements,  and 
can  be  applied  whenever  two  integers  can  be  stored  in  each  element  of 
the  control  array. 

The  iteration  is  now  carried  out  as  follows:  KONCO}  supplies  the  index  I 
of  the  first  part  to  be  considered  and  the  number  J  of  the  coefficient  set  to  be 
used.  We  then  compute  U(I)  as  follows: 

U(IJ  =  U(I)  +«|-U(I)  +  [ci(J)  ♦  U(I+kj)  +  C2(J)*U(I-kp  +  C3(J}*U(I+k2) 

+  C4(J)  ♦  UCl-k^)  +  C5(J)  ♦  UCI+k^)  +  C6( J)  ♦  UCl-k^) 

We  record  the  "error"  term,  which  is  all  of  the  right  hand  side  of  the  equation  except 
the  U(I)  term. 
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When  this  calculation  Is  completed,  we  pick  up  KC^(I)  which  gives  us  the 
next  values  of  I  and  J  .  When  KQN(I)  Is  zero,  this  Iteration  Is  complete, 
and  we  check  for  convergence  and  repeat  the  entire  Iteration  procedure  If 
necessary. 

The  entire  programming  procedure  Is  shown  In  detail  In  the  attached  block 
diagram  (Diagram  4.1).  The  diagram  should  be  self-explanatory,  but  notes  are 
appended  In  certain  places  to  aid  In  tracing  the  program  steps. 

In  order  to  abbreviate  the  diagram,  we  will  use  the  following  conventions: 

Symbol  "1“  will  Indicate  that  the  procedure  Is  to  be  carried  out  In  each  of 
the  six  directions.  That  Is,  the  statement  "Compute  means  compute  B^, 

B^,  Bj,  B^,  Bg,  and  B^  for  the  given  values  of  x,  y,  and  z  , 
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W  (U.A.) 
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5.  Examnlaa.  TO  support  our  contention  that  Method  D  is  superior,  from  the 
point  of  view  of  both  speed  and  accuracy,  to  any  other  method,  whether  discrete 
or  continuous,  typical  numerical  evidence  will  now  be  presented  by  means  of 
several  examples.  Each  example  was  run  on  the  CDC  1604  at  the  University 
of  Wisconsin. 

Example  1.  Let  G  be  the  spherical  sector  in  the  first  octant  bounded  by  the 

2  2  2 

surfaces  whose  equations  are  xsO,  ysO,  zsO,  x  -i-y  -t-z  >1  .  Let 

-  -  -  2  2 
F  s  0,  (x,  y,  z)  =  (0,  0,  0),  h  =  .1  and  4  s  x  +  2y  -  z  •  \  contains  410 

points.  Using  overrelaxation  with  a  zero  Initial  vector.  Method  D  yielded  results 

2  2 

correct  to  at  least  nine  decimal  places  to  the  exact  solution  u  =  x  4  2y  -  z  . 

The  running  time  was  1  minute  30  seconds. 

3  3 

Example  2.  Example  1  was  modified  by  setting  ^sxy-xy  -5z.  Selected,  but 

typical,  results  are  recorded  under  u^  in  Table  1.  The  running  time  was  1  minute 

3  3 

53  seconds  and  the  exact  solution  was  u^xy-xy  -  Sz. 

Example  3.  Example  2  was  modified  by  refining  the  grid  to  h  =  .  05  .  contained 
3721  points.  Selected,  but  typical,  results  are  recorded  under  u^^^  in  Table  1. 

The  running  time  was  8  minutes  11  seconds. 

Example  4.  Let  G  be  the  spherical  sector  in  Example  1.  Let  F  b  (see  [5]}  , 
*****  2 

(x,  y,  z)  =  (0,  0,  0),  h  =  .1  and  42x4  2y4z  .  contains  410  points. 

Method  D  was  applied  with  nonlinear  overrelaxation  [14]  and  with  zero  initial 
vector.  The  rurmlng  time  was  2  minutes  33  seconds.  Selected,  but  typical, 
results  are  recorded  under  u^^^  of  Table  2. 
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Examole  5.  An  extended  liebman  method  [10]  was  applied  in  Example  4  in  place 
of  nonlinear  overrelaxation.  Newton's  method  was  applied  to  solve  each 
equation.  The  running  time  was  4  minutes  20  seconds.  The  results  agreed  with 
those  of  Example  4  to  at  least  eight  decimal  places.  (With  regard  to  this  example, 
one  should  also  consult  [15]. } 

Example  6.  Let  G  be  the  four  dimensional  spherical  sector  defined  by 

2  2  2  2 

G  =  X  >  0,  y  >  0,  z  >  0,  «  >  0,  X  +y  +z  +«  <1}  .  Let 

2222 

F  a  0,  (x,  y,  z,  <i) }  s  (0,  0,  0,0),  h  =  .  1  and  #*x  -2y  +3z  -Zu  .  contained 
803  points.  The  difference  analogue  of  the  four  dimensional  Laplace  equation, 
that  is,  the  extension  of  (2.1)  to  four  dimensions  can  be  written  for  this  problem 
in  the  form 


(5.1) 


4 

Z 

i=l 


-u. 


2i-l 


IL 


^2i-l*'2i  *'21-1^^2i-l'‘’*'2i^  ^2i^^21-l'*'**2p 


s  0 


Using  overrelaxation  with  initial  vector  the  zero  vector,  Method  D  approximated 

2  2  2  2 

the  exact  solution  u  =  x  -  2y  +  3z  -  2(i)  at  each  point  to  at  least  nine  decimal 
places.  The  running  time  was  3  minutes  58  seconds. 


Example  7 .  In  this  example  we  will  show  how  Method  D  can  be  applied  even 
when  G  is  not  simply  connected.  Let  Sj  be  the  "outer"  oubic  surface  with 
vertices  (0,  0,  0),  (2,  0,  0),  (0,  2,  0),  (0,  0,  2),  (2,  2,  0),  (2,  0,  2),  (0,  2,  2), 

(2,  2,  2)  and  let  S_  be  the  "inner"  cubic  surface  with  vertices  (j,  ^),  (7,7,7  ) , 

M  fc  b  to 


'2'  2’  2  v2>  2  *  2^'^2’  Z*V*^Z*  Z*  2”  2’  2* ’'2*  Z*  2  * 


as  shown  in  Diagram  &  1. 
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Let  R  be  the  region  between  and  •  Let  and  . 

Also,  it  is  given  that  1  on  S^,  #a  o  on  and  paQ  .  Then,  setting 
(x,  y,  z)  =  (0,  0,  0)  and  h  =  .1,  Method  D  was  applied  to  the  resulting  Dirichlet 
problem.  contained  784  points.  The  symmetry  of  the  solution  was  incorporated 
into  the  program  so  that  only  one  eighth  of  G  had  to  be  considered.  Over- 
relaxation  was  applied  with  initial  vector  the  zero  vector  and  selected,  but  typical, 
results  are  recorded  under  u^^^  of  Table  3.  The  running  time  was  2  minutes 
26  seconds. 

Example  8.  Example  7  was  modified  by  refining  the  grid  to  h  =  .  05  .  R^^  contained 
6669  points.  Selected,  but  typical,  results  are  recorded  under  u^^^  in  Table  3. 

The  running  time  was  12  minutes  49  seconds. 
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TABLE  1 


Approx.  Sol.: 
ut‘> 

Approx.  Sol.: 

(2) 
u'  ' 

Exact  Sol. : 

3  3,  , 

u=x  y-xy  +5z 

-0.  50000000 

raHii 

-.  50000000 

-3. 50000000 

-3. 50000000 

-3. 50000000 

-1. 00059993 

-1.00060000 

-1.00060000 

-0. 50239984 

-0. 50239998 

-0. 50240000 

-3. 00239930 

-3.00239992 

-3. 00240000 

-1. 00599922 

-1.00599990 

-1.00600000 

-4. 00599842 

-4. 00599974 

-4. 00600000 

-1.02099693 

-1.02099962 

-1.02 100000 

-0. 53359660 

-0. 53359961 

-0. 53360000 

-3.5336 1012 

-3.53360023 

-3. 53360000 

-2.05039221 

-2.05039855 

-2. 05040000 

-0. 49940004 

-0.49940000 

-0.49940000 

-2.99940018 

-2.99940002 

-2.99940000 

-0. 50000000 

-0. 50000000 

-0. 50000000 

-3.00000000 

-3. 00000000 

-0.50299980 

-0. 50299997 

-0. 50300000 

-0.50959933 

-0.50959992 

-0. 50960000 

-3.00959709 

-3.00959968 

-3. 00960000 

-2. 52099448 

-2.52099937 

-2.52100000 

-1.03839447 

-1.03839930 

-1.03840000 

-3.03838782 

-3.03839868 

-3.03840000 

-3. 06297700 

-3.06299759 

-3.06300000 

-0.497600 16 

-0.49760002 

-0.49760000 

-2.99760070 

-2.99760008 

-2.99760000 

-0.49700020 

-0.49700002 

-0.49700000 

-2,99700092 

-2.997000 10 

-2.99700000 

-0.50000000 

-0. 50000000 

-0. 50000000 
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TABLE  1  (Continued) 


Approx.  Sol.: 
u“> 

Approx.  Sol.: 

(2) 

u' 

Exact  Sol. : 

3  3  , 

u=x  y-xy  +5z 

-4. 50000000 

-4. 50000000 

-4. 50000000 

-4.008  39981 

-4.00839951 

-4. 00840000 

-3.02399285 

-3.  02399923 

-3.02400000 

-2. 04858917 

-2. 04859867 

-2.04860000 

-1.58398473 

-1.58399789 

-1.58400000 

-1.13198012 

-1.13199693 

-1.  13200000 

-1.19437299 

-1.  19439290 

-1.  19440000 

-1.99400130 

-1.994000 16 

-1.99400000 

-3.99400158 

-3.99400026 

-3.99400000 

-3.99040373 

-3.99040039 

-3.99040000 

-2. 49 1602  17 

-2.49160024 

-2.49160000 

-0.  50000000 

-0. 50000000 

-0.50000000 

-1.01799753 

-1.01799967 

-1.01800000 

-2.04798864 

-2. 04799869 

-2.04800000 

-1.59238087 

-1.59239767 

-1.59240000 

-2.  15360595 

-2.  15359356 

-2. 15360000 

-3.48800346 

-3.48800042 

-3.48800000 

-2.47900552 

-2.47900062 

-2.47900000 

- 1. 47600437 

-1.47600055 

-1.47600000 

-1.48200338 

-1.48200042 

-1.48200000 

-0. 50000000 

-0. 50000000 

-0. 50000000 

-3. 50000000 

-3. 50000000 

-3.50000000 

-1.08398798 

-1.08399811 

-1.08400000 

-0,97900307 

-0.97900037 

-0.97900000 

-2.97900485 

-2.97900064 

-2.97900000 

-1.96160854 

-1.96160104 

-1.96160000 

-0.95140680 

-0.95140090 

-0.95140000 
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X 

y 

z 

Approx.  Sol.: 
u“' 

Approx.  Sol.: 

(2) 
u'  ' 

Exact  Sol. : 

3  3  . 

u=x  y-xy  +52 

.6 

.3 

.6 

-2.95141421 

-2.95140156 

-2.95140000 

.6 

.4 

.5 

-2.45201490 

-2. 45200 137 

-2.45200000 

.6 

.5 

.5 

-2.46700303 

-2.46700050 

-2.46700000 

.6 

.6 

.4 

.7 

.  1 

.4 

-1.96640694 

-1.96640094 

-1.96640000 

.7 

.3 

-1.43701210 

-1.43700 156 

-1.43700000 

.7 

.2 

-0.91601231 

-0.91600 171 

-0.91600000 

.7 

.  1 

-0.40760854 

-0. 40760 124 

-0. 40760000 

.7 

-2.40764318 

-2. 407603 12 

-2.40760000 

.7 

.5 

-2.41599728 

-2.41599952 

-2.41600000 

.7 

.  1 

-0. 50000000 

.8 

.  1 

-0.  40401347 

-0.40400 140 

-0. 40400000 

.8 

,4 

-1.86803535 

-1. 86800377 

-1.86800000 

.8 

.  1 

-0. 34404041 

-0. 34400209 

-0. 34400000 

.9 

.  1 

.4 

-1.92799112 

-1.92800 167 

-1.92800000 

.9 

.4 

.  1 

-0. 26598085 

-0.26599834 

-0. 26600000 
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TABLE  2 


B 

D 

H 

u“> 

fl 

B 

z 

ut‘> 

.  1 

.  1 

.  1 

. 33379803 

.4 

.  1 

.77843093 

.  1 

.2 

.4 

.73417575 

.4 

.2 

.  1 

. 88700006 

.  1 

.3 

.5 

1. 04443022 

.4 

.3 

.3 

1. 28020294 

.  1 

.4 

.3 

1.08918215 

.4 

.3 

.7 

1. 6 1368477 

.  1 

.5 

.2 

1.23012008 

.4 

.4 

.4 

1. 58080909 

.  1 

.6 

.4 

1.57000816 

.4 

.5 

.7 

1.95821289 

.  1 

.7 

.4 

1.75731038 

.4 

.6 

.3 

1. 89  143826 

.  1 

.8 

.2 

1.81786813 

.4 

.7 

.5 

2.  12651462 

.  1 

.9 

.4 

2.06963333 

.4 

.8 

.2 

2. 13369786 

.2 

.  1 

.5 

.72353422 

.5 

.  1 

.4 

.95402334 

.2 

.3 

. 8020349 1 

.5 

B 

.3 

1.  14019526 

.2 

.9 

1.44682719 

.5 

D 

.2 

1.29857485 

.  2 

.9 

1.63582348 

.5 

.4 

1.66722297 

.2 

.8 

1.71152454 

.5 

B 

.7 

1.99755020 

.2 

.5 

.8 

1.87751902 

.5 

H 

.2 

2. 06 195601 

.2 

.7 

.  1 

1.70257965 

B 

.  1 

.3 

.97486870 

.2 

.8 

.  1 

1. 88993511 

B 

B 

.6 

1.45675948 

.3 

.  1 

.  1 

.55309066 

B 

H 

.3 

1.66927338 

.3 

.  1 

.6 

.94374260 

.6 

.6 

1.98270586 

.3 

.2 

.7 

1.  30608 156 

.6 

.3 

2.  13505733 

.3 

.3 

.5 

1.33513052 

B 

fl 

.7 

1.39334550 

.3 

.4 

.4 

1.47053424 

H 

.2 

1.5957  1747 

.3 

•5 

.2 

1.  5 16 16497 

fl 

fl 

.3 

1.89236582 

.3 

.6 

.  3 

1.79  105773 

.8 

.  1 

.2 

1.09268284 

.3 

.6 

1.85851258 

.8 

fl 

.5 

1.48  122715 

.3 

.7 

.4 

2.01526321 

.8 

fl 

1.75110486 

.3 

.9 

.  1 

...M 

2.1489  1247 

.  1 

.3 

1.21296598 
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TABLE  3 


X 

y 

z 

(2) 
u'  ' 

1.0 

1.0 

1.6 

.21318413 

.21340617 

1.0 

1.4 

1.6 

.25474801 

.26000151 

1.0 

1.6 

1.3 

. 23039266 

.23203625 

1.0 

1.8 

1.9 

.93922824 

. 94029832 

1.0 

1.9 

1.6 

. 87962760 

.88193903 

1. 1 

1.5 

1.7 

. 52860578 

.  53955766 

1. 1 

1.6 

1.6 

.467302 12 

. 477967 16 

1.1 

1.8 

1.5 

.70210223 

.  708 15556 

1.1 

1.9 

1.0 

.  8 1346 190 

.81395052 

1.2 

1.  1 

1.9 

.  8 18 12623 

.  8 1902279 

1.2 

1.6 

1.  1 

. 22076057 

.22 147398 

1.2 

1.7 

1.  1 

.43191770 

.  43329 169 

1.2 

1.8 

1.  1 

.63036172 

.63180745 

1.2 

1.9 

1.7 

.91165183 

.91351895 

1.3 

1.5 

1.7 

.54872314 

.  56 194398 

1.3 

1.7 

1.5 

.54872314 

.  56  194398 

1.3 

1.8 

1.2 

.64900339 

.65215013 

1.4 

1,0 

1.7 

.47342317 

.479  33409 

1.4 

1.2 

1.8 

.67122668 

.67620583 

1.4 

1.6 

1.4 

. 29590727 

.  307046 19 

1.4 

1.8 

1.9 

.94863874 

.95008233 

1.5 

1.2 

1.6 

.31897  168 

.33873950 

1.5 

1.6 

1.5 

.41057277 

.45137840 

1.5 

1.8 

1.4 

.73764103 

.74655664 

1.5 

1.9 

1.6 

.90921585 

.91255196 

1.6 

1.  1 

1.8 

.75512719 

.76003480 

1.6 

1.2 

1.9 

.883187  14 

.88577147 

TABLE  3  (Continued) 


X 

y 

z 

u“> 

(2) 
u'  ' 

1.6 

1.5 

1.6 

.56670759 

. 58864922 

1.6 

1.7 

1.6 

.7599  1236 

. 76967969 

1.6 

1.9 

1.3 

.88834726 

. 89 126535 

1.7 

1.0 

1.0 

.42086795 

. 42 13485 1 

1.7 

1.2 

1.0 

.42989416 

.43111151 

1.7 

1.4 

1.3 

. 496 16024 

.50444129 

1.7 

1.6 

1.5 

.7008 1037 

.71379555 

1.7 

1.8 

1.5 

.86041861 

.86533931 

1.8 

1.  1 

1.3 

.64295832 

.64555233 

1.8 

1.2 

1.7 

.  8202 1476 

,82397741 

1.8 

1.4 

1.4 

.70371313 

. 7  1 143949 

1.8 

1.6 

1.7 

.88560303 

. 8897  167  1 

1.8 

1.9 

1.3 

.94416441 

.94549284 

1.9 

1.0 

1.4 

.83606017 

.83822357 

1.9 

1.3 

1.  1 

.82513613 

.8266 1512 

1.9 

1.3 

1.7 

.91578506 

.91787661 

1.9 

1.5 

1.7 

.93179483 

.93410454 

1.9 

1.7 

1.0 

.90873890 

.91041024 

1.9 

1.9 

1.0 

.96974268 

.97026699 
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